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A lattice calculation of the branching ratio for B K*j 

UKQCD collaboration, presented by Hugh Shanahan ^ 

^Department of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3JZ, Scotland 

Preliminary results are presented on the calculation of the relevant form factor for the radiative decay B K*j . 
We find that the form factor is weakly dependent on the spectator quark mass. We compare the results obtained 
for a full chiral extrapolation and where this independence is assumed. Both results are found to be statistically 
consistent with the experimental data. 
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1. Introduction 

The decay B K*j is unique in being simul- 
taneously of interest to experimentalists, and the 
of supersymmetry and lattice field theory commu- 
nities. As the quark decay h ^ sj cannot occur at 
tree-level, it is sensitive to a range of parameters 
in the Standard Model, including iVtsT^jj and mj. 
It is also sensitive to new physics from two Higgs' 
doublets models [1], supersymmetry [2], and even 
supergravity [3]. While there is currently debate 
on how much could be determined from this de- 
cay for these models, it is clear that h ^ sj is 
a good probe of possible new physics, and that 
accurate theoretical predictions are essential. 

B K*j is an excellent candidate to ex- 
plore this quark decay experimentally, as it has 
a clean signature. The CLEO collaboration have 
published a branching ratio for B K*j of 
(4.5 ± 1.5 ± 0.9) X 10-5[4], albeit with quite 
low statistics. However, unlike h ^ sj , there 
are large soft-gluonic contributions to this de- 
cay, and different models have given a wide range 
of different answers. It is clear then that non- 
perturbative methods, and lattice field theory in 
particular offer a way ahead. 

Following the notation of Grinstein, Springer 
and Wise [5], the branching ratio can be evaluated 
in the leading log approximation from the matrix 
element 

'-^0Crim,)VuV,U^*{K*\sa,^q-'bn\B) 

where C'rlmi,) is calculated perturbatively and 
qf' , rj^ are respectively the momentum and polar- 
isation of the photon. Bernard, Hsieh and Soni 



[6] demonstrated the feasibility of a lattice cal- 
culation of [B\h<T fii^q" s\K*) , which can expressed 
using three form factors Ti^q"^) multiplying three 
possible Lorentz structures C^: 



(B\ba,^q^s\K*) = J2c;T,iq'). 



(1) 



For on-shell photons, q'^=0 and the contribution 
from the third form factor vanishes, furthermore 



T2(g2=0) = -iTi(g2=0). 



(2) 



Hence, the branching ratio can be expressed as 
a product of a number of parameters which are 
quite accurately known, C7(mj)^, which depends 
on mt, and Ti(q'^=0)^ . This paper concentrates 
on evaluating Ti(q'^=0). 

2. Computational details 

The analysis for this calculation is based on 60 
quenched configurations for a 24^ x 48 lattice at 
l3 = 6.2, generated on the Meiko i860 Computing 
Surface based at Edinburgh. From string tension 
measurements [7], the lattice spacing a was deter- 
mined to be 2.73(5) GeV. The algorithms used in 
the gauge field generation and calculating matrix 
elements of the form {A\qiTq2\B) , where \A) and 
\B) are meson states, are explained in [7] and [8]. 

The quark propagators are evaluated using 
the 0(a) improved ("clover") fermionic action, 
with periodic boundary conditions. To simulate 
the spectator (light) quarks, the propagators are 
evaluated at ki = 0.14144,0.14226 and 0.14262, 
whose pole masses are approximately that of the 
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Figure 1. Ti^q"^), with a linear fit. Tlie dotted 
lines represent the 68% confidence levels of the 
fit at q^=0. 

strange quark, given that kP'^^" = 0.1419(1) [9]. 
For the other quark in the K* state, the first two 
K values are used as they straddle the physical 
strange quark mass. Heavy quark propagators 
are computed at Kh = 0.f29 and 0.f21, with pole 
masses of approximately 1.9 and 2.5 GeV. The 
operator haf^^s is improved to 0(a). 

Gauge invariant smearing is used for the ex- 
tended pseudoscalar field to improve the signal. 
We apply time-reversal symmetry in the evalu- 
ation of '^^£p{V(k,£)\'safivh\P(p)) for the same 
reason. The matching of the lattice operator with 
the continuum [10] has a negligible contribution 
and is not considered here. 

Ti is evaluated from the above by considering 
different components of the relation 

4(A;V -p«A;'')Ti(g2) = 

s"PP'^J2e,{V(k,e)\sa,,h\P(p))r/. (3) 

The resulting summed data is averaged on both 
sides of the lattice. A correlated fit to a constant 
is carried out for timeslices 11 to 13. The statis- 
tical analysis for this and all the subsequent fits 
of this paper is based on the bootstrap algorithm, 
using 250 subsamples. 

The form factor is evaluated for seventeen dif- 
ferent values for the spatial momentum q, rang- 
ing in magnitude from to 27r/(12a), and for 
p = (7r/(12a),0,0) and (0,0,0). For this choice 
of momenta and k values, the range of q"^ always 
includes zero and is small in comparison to the 



scale expected from the vector meson dominance. 
Hence, for each combination of and Kh, 

Ti^kijksjhhiq"^) is linearly interpolated to q'^=0. 
An example of this fit is shown in figure 1. 

The variation of Ti with respect to the quark 
masses is now considered. 

3. Extrapolation to physical mass scales 

The data set Ti(ki, Kg, Kh',q'^=0) is used to 
derive the form factor at the appropriate mass 
scales. Given the high statistics of this calcula- 
tion, Ti(ki, Kg, Kh',q'^=0) can be extrapolated to 
the chiral limit. It is fitted to the function 

Tl (ki ,Ks,Kh; q^ = 0) = Tl (Kcrit ,Ks,Kh; q^ = 0) 

+ A(Ks,Kh)mg light, (4) 

where niqUght is the pole mass of the light quark. 
An example of such a fit is in figure 2. The re- 
sult of this fit, Ti(Kcrit, Kg, Kh',q'^=0), is linearly 
interpolated to KPj^y^ . The two remaining values 
of the form factor are extrapolated to the B mass 
scale with the function 

Ti(«:er,t,«:ft;g'=0) = Tf''*- + ^, (5) 

giving a result of Ti{q'^=Q) = Q.l'ot^^^- We note 
that the chiral extrapolation of the form factor ex- 
hibits a very weak dependence on the light quark 
pole mass. If we assume the data is indeed in- 
dependent of the light quark mass and fit to a 
constant, we find the per degree of freedom to 
be of the same order as fitting the data to a linear 
function. If the data from this fit is continued to 
the physical regimes one finds Ti((;^=0) = 0.15^4, 
significantly reducing the statistical errors and 
maintaining the central values. Both extrapola- 
tions to the B mass scale are shown in figure 3. 

For comparison, the branching ratio of 5 ^ 
A'* 7 from CLEO is converted into an experi- 
mental value for the form factor, Us- 
ing the parameters from the particle data book, 
and taking mj = 4.39 GeV, T^'^p^ is found to be 
0.23(6), 0.21(5) and 0.19(5) for rut = 100, 150 
and 200 GeV. These results are statistically con- 
sistent with those found from our analysis. 



3 



T,, chiral limit 



/Cj,=0. 13900 
« =0,14144 



0.0 ^ 
0.00 



0.08 0.03 

"Iq light 



Figure 2. Chiral extrapolation of Ti(q'^=0). 
rriq light is the lattice pole mass. 



4. Conclusions 

The calculation of matrix elements from three 
point functions on the lattice is rapidly approach- 
ing maturity. In our case, the derivation of the 
matrix element {V\s<t fj_uh\P) in the quenched ap- 
proximation presents no conceptual difficulties 
and gives a good signal. The question remains 
whether extrapolating this data to the chiral limit 
and the appropriate heavy mass scale is valid, as 
the former is susceptible to large systematic er- 
rors due to finite volume effects and the latter 
does not have a clear theoretical justification for 
the ansatz used. 

In the first case, it is very encouraging that the 
form factor Ti has no, or at least a very weak 
dependence on the light quark mass, indicating 
a stability in the data. To address the second 
issue, we are currently repeating the calculation 
on the same lattice, using two new heavy quark 
masses at = 0.125 and 0.133. Furthermore, 
we plan to repeat the calculation on a 16^ x 48 
lattice at /3 = 6.0, generating the heavy quark 
propagator with a hopping parameter expansion. 
This will allow us to explore a large number of 
heavy quark masses. It shall also check for any 
systematic errors due to the lattice spacing. 

With our current data, we shall repeat the anal- 
ysis presented by Soni at this conference [11], 
where the form factor T2 is extracted. This can be 
used as a further determination of the systematic 
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Figure 3. Extrapolation of Ti((;^=0) to m^. The 
fitted points are displaced for clarity. 
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